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Description 

Field of the invention. 

s Th present inv ntion is in the field of digital radiography. The invention more specifically relates to a method of 

recognising an irradiation field in a system wherein a radiation image is recorded on a photostimulabte phosphor screen 
by means of an irradiation field stop (i.e. means for partially shielding against irradiation). 

Description of prior art 

10 

Radiologists may protect their subjects against unnecessary exposure to X-rays by use of X-ray opaque ('collima- 
tion') material. The material is placed in the path of the X-ray beam so as to shield those areas of the patient which 
are not regarded as diagnosticalry important. Besides reducing patient dose this technique has the additional advantage 
of limiting the amount of scattered radiation in the image field of interest. The regions of the resulting image which 
is comprise the shadow cast by the X-ray opaque material (shadow regions, cfr. example in figure 1) are exposed only 
by scattered radiation. 

In computer-assisted radiology, the X-ray image is digitised as a two dimensional array of numbers, whose mag- 
nitude is related to the intensity of X-rays arriving at the detector, via the patient. The values in the array may be 
reseated, prior to display or printing, in order to maximise the visual contrast in the area of diagnostic interest. Such 
20 rescaling depends, in part, on the intensity histogram of the X-ray image. To achieve this function automatically, best 
results are obtained if the 'shadow* cast by the collimation material is excluded from consideration, since it cannot 
contain any useful information, but would otherwise dominate the tower end of the intensity histogram. 

The presence of the collimation shadow region can cause an additional problem in the display of the image on 
film. Therein, the shadow region is relatively bright, and if displayed unmodified, may impair diagnosis of subtle lesions 
25 due to dazzle, especially if the unexposed region is relatively large. 

The U.S. patent 4,952,807 by Adachi addresses the first problem, wherein a method is disclosed tor adjusting the 
image processing conditions by selecting only those pixels which correspond to the object or to the background, i.e. 
the 'signal' region, under the assumption that the collimation material does not partly cover the object, i.e. the shadow 
region does not touch the object of the image, but rather its surrounding background. In many kinds of examinations 
30 however, parts of the patient will be shielded from exposure to avoid useless X-ray load, as is the case with lumbar 
spine and shoulder examinations. 

The U.S. patent 4,804,842 from Nakajima also addresses the first problem, wherein a method is disclosed for 
removing the lower part of the histogram corresponding to the pixels in the shadow region in order to adjust the image 
processing conditions. 

35 However, techniques for histogram modification alone, do not address the second problem, i.e. that of dazzle. Its 

solution demands an explicit detection of the shadow region, as distinct from the signal region. The offending shadow 
region may then be excluded from the display process. In addition, the first problem of rescaling for optimum visual 
contrast can then proceed on the signal region alone, without the complicating influence of the shadow region on the 
intensity histogram. 

40 Previous work in the field of automatically delineating the signal/shadow regions has relied on piecemeal attempts 

at the solution. 

Many of the proposed methods are based on collecting candidate signal/shadow border pixels (i.e. edge pixels at 
the interface between signal and shadow regions) by thresholding on first differences along predetermined lines that 
might cross the signal/shadow border. A typical example of this approach is proposed by Funahashi in the U.S. patent 

45 4,970,393. The signal region is recognized as a region enclosed by pixels where the first difference exceeds some 
specified threshold, either positive or negative. 

A refinement of this technique is proposed by Takeo in the European patent 0 342 379, where multiple candidate 
pixels along a searching line are ranked, and only the highest rank candidates on each searching line are considered 
in assembling the signal/shadow boundary. In both methods the underlying assumption is that the signal/shadow bound- 

so ary has strong contrast everywhere. In many cases however, contrast aiong parts of the signal/shadow boundary is 
nearly zero, especially if the collimation material is partially covering very dense parts of the body. So if the recognized 
signal/shadow boundary is assembled from high difference pixels many gaps will occur in general. 

In the European patent 0 285 174 by Funahashi a straightforward application of the basic Hough transform tech- 
nique is proposed to delineate the signal/shadow boundary, which is assumed to be polygonal. The Hough transform 

55 is applied to prospective edge points obtained by differentiation. 

The underlying assumptions are that th signal/shadow boundary is piecewise linear, and has high contrast with respect 
to the contrast of line features within the signal region. However, experimental evidence shows that edges within the 
signal region, e.g. bone edges, frequently have higher contrast than the signal/shadow boundary. As a consequence 
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points in Hough transform space that correspond to straight bone boundaries may have higher scores than some of 
the actual signal/shadow borders, in which case the subsequent delineation procedur will fail. 

A diff rent approach is disclosed in the European pat nt 0 360 231 by Funahashi. The image is subdivided into 
mutually adjacent non-ov flapping tiles. A statistical parameter is computed indicating the edge content of the tile, 

s either by differentiation in image space, or in frequency space. The tiles are accordingly classifi d as b longing to the 
relevant image region or not. 

We ran elaborate experiments and concluded that single criteria based on edge contrast, or local statistical pa- 
rameters at both sides of a prospective signal/shadow boundary, such as mean value, min, max, variance or spatial 
frequency coefficients are inadequate in delineating the irradiated field with a high level of confidence, say > 98 percent. 

10 We found that the signal/shadow border contrast will vanish in some cases; anatomical structures may be straight and 
have high contrast edges; the image signal (which represents local dose) in the shadow region will occasionaly be 
higher than the signal in the unshielded (i.e. diagnostic) image region - in cases where the scattered radiation is high 
and some parts of the body are difficult to penetrate by radiation. Also we found that there is no essential distinction 
between image areas corresponding to very poorly penetrated bone, and areas corresponding to shielded parts. Both 

is kinds of areas will be characterised by an extremely low mean value (i.e. dose), similar noise level (both quantum 
limited), and by the absence of other features (since the signal to noise ratio is very low). 

In order to obtain a high success rate, a preferred method should consider many hypotheses as to the location of 
the signal/shadow boundary, and keep or reject each hypothesis 'as a whole", instead of keeping or rejecting individual 
components of a hypothesis, such as edges. 

20 Multiple hypotheses as to the location of the signal/shadow boundary are considered in the method disclosed by 

Tanaka in the U.S. patent 4,952,805. For each of the corresponding candidate signal regions, two signal histograms 
are computed: one for the pixels inside the candidate irradiation field and one for the pixels outside. The candidate 
boundary which is characterised by the largest degree of interclass separation as determined from both histograms is 
selected to be the recognised irradiation field. 

25 The degree of interclass separation is defined in such a way that a high value will result if there is very little overlap 
between both histograms. Hence this criterion will be effective only if the range of signal values within irradiation field 
are well separated from the signal values in the shielded region, which does not hold in many practical cases as we 
verified on the basis of digital radiographs from various examinations. 

Also, there is a practical limitation regarding the usability of the above mentioned approach, in that the number of 

30 hypotheses as to the location of the irradiation field may be prohibitively large. 

Previous work in general has focussed on methods in which finding the boundary of an irradiation field essentially 
relies on a very small number of mostly local criteria. Reviewing a wide variety of examinations confirms that individual 
criteria will only work in a very limited number of cases. A high success rate can only be achieved if many sources of 
evidence are combined, ranging from local image features to global constraints. 

35 The behaviour of the radiologist, in placing the collimation material around the patient, may broadly be characterised 

by a set of simple rules. This invention addresses only the use of collimation material which casts a shadow whose 
outline is piecewise linear. More specifically, the boundary between signal, and shadow regions is assumed to be a 
rectangle of arbitrary location and orientation. It is assumed that the signal region is on the interior, and that any part 
of the rectangular boundary may, or may not, have a projection onto the (rectangular) detector, and therefore the image 

40 array; see figure 1. 

Objects of the invention 

It is an object of the present invention to provide a method of recognising an irradiation field in a system wherein 
45 a radiation image has been recorded by using an irradiation field stop (collimation material). 

It is a further object to provide a method for determination of the location of the signal/shadow boundary of an X- 
ray image, wherein the x-ray image is represented by a digital signal representation and wherein a signal/shadow 
boundary is the locus of interface between the signal and the shadow regions resulting from the obscuring by the 
collimation material. 

50 it is a still further object to provide such a method in a system wherein a radiation image is recorded in a photo- 

stimulable phosphor screen and wherein said radiation image is read by scanning the screen with stimulating radiation 
and by detecting the light emitted upon stimulation. 

it is stilt a further object of the present invention to aid the process of contrast-enhancement of an X-ray image, 
by permitting such to perform its calculations based solely on the diagnostically useful signal part of the X-ray image, 

55 and to reduce the brightness in shadow region to avoid dazzle during diagnosis. 
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Statement of the Invention 

The present invention provides a m thod for the determination of the location of th signal/shadow boundary in 
an X-ray image represented in a digital signal representation comprising the st ps of : 

5 

i) Exacting low-level primitives from the X-ray image {X(i,j)} t ii) Forming a reduced number of intermediate-level 
primitives from the low-level primitives, 

iii) Building hypotheses being possible solutions regarding the location of the signal-shadow boundary from com- 
binations of intermediate-level primitives, during which each combination is subject to intermediate-level tests, 

10 whereupon partial or complete hypotheses are rejected or accepted, 

iv) Performing high-level verification tests on each hypothesis, whereupon hypotheses are rejected, or accepted 
at some cost said cost defining a ranking of the acceptability of a hypothesis, 

v) Selecting the hypothesis with the least cost. 

1$ The method of the present invention is a major departure from earlier methods in that, in general, it generates 

many hypotheses as to the location of the boundary. The final decision as to the correct hypothesis is contingent on 
the results of a number of high-level tests applied to each. 

The method permits codification of expert-system type knowledge as to the geometry, topology, and statistics of 
the signal/shadow boundary, and associated regions. 

20 The invention is, for example, applicable to a digital image signal obtained by exposing a photostimulable phosphor 

screen to an image originating from a partially shielded exposure of an object to X-rays, by scanning the screen with 
stimulating irradiation, detecting the light emitted upon stimulation and converting the detected light into a digital signal 
representation. 

Commonly the digital signal is subjected to image processing such as contrast enhancement and signako-density 
2S mapping before being applied to an output device for reproduction or display The result of the process of contrast- 
enhancement of an X-ray image is enhanced by permitting such to perform its calculations based solely on the diag- 
nostically useful signal part of the X-ray image, and to reduce the brightness in shadow region to avoid dazzle during 
diagnosis. 

In the statement of the invention and in the description hereinafter the term "X-ray image ■ is used in the meaning 
30 of a two-dimensional digital image repesentation as an array of numbers, the magnitudes of which are directly related 
to the intensity of the X-rays arriving at a 2-dimensional physical detector such as a photo-stimulable phosphor screen. 

By "primitives" is meant features, such as edges, which are extracted from the image as the result of low-level 
image processing, cfr. fig. 2. 

An "intermediate-level primitive" is a combination of (basic) primitives which may be considered collectively when 
35 building a hypothesis. Intermediate-level primitives are generally fewer, and (spatially) larger than are (basic) primitives. 
Each primitive is constructed from a mutually exclusive set of pixels. Each intermediate-level primitive is constructed 
from a mutually exclusive set of (basic) primitives. 

By the term "signal region" is meant the spatial area of an image which is not obscured by collimation material 
from a direct line-of-sight of the X-ray source, see figure 1b. 

40 

The "shadow region" is the area of the image which is obscured by the collimation material from a direct line-of- 
sight of the X-ray source; see also figure 1b. 

The "signal-shadow boundary" is the locus of the interface between the signal and the shadow regions. E.g. the 
pixels at the interface between the signal and the shadow regions. 

45 

A low-level primitive can be a line - being a data structure to represent an edge, and associated data, including its 
parameterisation - or a tile - being a small rectangular area of pixel resulting from a tessellation of the image array. 

Steps i) and ii) of the method of the invention achieve a reduced set of primitives which can be assembled into 
hypotheses as to the location of the signal-region. 
50 a "hypothesis" is defined as a possible solution posited by the method of the present invention, said solution being 

a segmentation (division into similar regions) of the original image into signal and shadow regions, e.g. a two-dimen- 
sional array of numbers isomorphic with the X-ray image, with pixel labels for signal and shadow regions. 

Hypotheses are generated from the reduced (intermediate-level) set because, in general, the number of low-level 
primitives is too numerous to perform an exhaustive evaluation (step iv) of every possible combination (step iii). 
55 Specifically, if the low-level primitives are lines, then the intermediate-level primitives would be line-clusters, the 

result, for instance, of a linear agglomerative clustering process. 

A line cluster is then one or more lines which are deemed to be in some manner similar, and hav been accordingly 
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grouped together. 

E.g. one or mor lines which have been grouped together by a linear agglomerative clust ring process. In th 
example of figure 3, the lines 1b and 1c may be clustered togeth r to form a line clust r {1b,1c}. Similarly, in the 
same example, the likely clusters are {1a}, {1d}, {1e, 1f, 1g}. 

5 

Alternatively, a low-level primitive can be a tile, i.e. a small rectangular array of pixels resulting from a tessellation 
of the image array In this case an intermediate-level primitive would be a tile-cluster, being a group of spatially con- 
nected tiles having in common a feature or an attribute. 

Step iii) combines the intermediate-level primitives into complete hypotheses posited by the system. 

10 

The definition of the word "hypothesis" was given hereinbefore, by which term is meant a possible solution regarding 
the shape and position of signal and shadow regions. 

Such a solution is a segmentation or division of the image into signal and shadow regions. 
One of the possible representations of such a segmentation is a two dimensional array of numbers, isomorphic 
is with the X-ray image, the values of which represent a unique label which identifies the corresponding region in the 

X-ray image as either signal or shadow region. 

A label generally is a pixel-value which associates the pixel, i.e. the location in the image array, with a particular 
class (having a particular quality or attribute). 

20 in order to limit the number of hypotheses to be submitted for evaluation (step iv), the method includes provision 

for tests which may be applied during the construction of the hypothesis, which serve to detect and reject, at an early 
stage, an incorrect hypothesis. 

Specifically, if the intermediate-level primitives are line-clusters for instance, then a possible test is on the mutual 
orthogonality of the line-clusters as putative sides of the signal/shadow border. Thus combinations of line-clusters 

25 which are not going to form a rectangle can be rejected at an early stage. 

Step iv) permits the evaluation of hypotheses successfully assembled from the intermediate-level primitives. A 
major feature of this invention is the generation of multiple hypotheses, where in each hypothesis is a candidate for 
the complete signal/shadow boundary. The high-level tests permit evaluation of each hypothesis, specifically of those 
attributes which can only be determined once a complete signal/shadow boundary is hypothesised. By such high-level 

30 tests, hypotheses at this level may be ranked (or possibly rejected) in order that a single candidate may eventually be 
chosen. 

Specifically, the ranking of hypotheses may be based on a scalar quantity representative of either the net positive 
or negative evidence in favour of, or against that hypothesis. The distinction is immaterial: the highest ranking hypoth- 
esis is that most like the idealised model, in the former case this will have the largest (positive score), in the latter: the 
35 smallest negative score. 

The evidence may be based on any of a number or combination of primitives, and less local features of the original 
X-ray image. For instance, negative evidence may be accumulated against a candidate signal/shadow boundary for 
each corner of the signal/shadow boundary according to (a monotonic function of) the magnitude of the difference of 
the angle at the comer from 90 degrees. Such a test would attempt to enforce the rectangular geometry of the model. 
40 Broadly, tests may seek to enforce the geometry, topology region statistics, and edge statistics required by the model. 

Within these categories, other possible sources of evidence may involve tests for: 

Edge statistics: 

45 - Linearity: The straightness for each of the candidate signal/shadow boundary sides (however they were achieved). 

Existence: The accumulated evidence along the length of each side, computed from an edge-image, such as that 
furnished by Canny edge-detection for instance. 

Continuity: The length of each candidate signal/shadow boundary side for which there exists no evidence from the 
edge-detection process. 

so - Explanation: The ratio of used evidence to the total evidence. In this context, the evidence is the number of edge- 
pixels. 

The used evidence is the number of edge pixels in the signal/shadow boundary. 

The unused evidence is the number of edge pixels which lie on edges outside of, but colinear with, those edges deem- 
55 ded to be in the candidate signal/shadow boundary. 
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Geometry: 

Proportion: The aspect ratio of the rectangle (perhaps favouring rectangles with an aspect ratio similar to that of 
the X-ray image). 

s - Area: The area of the candidate signal-region (perhaps penalising very small rectangles). 

Position: The placement of the centre of the candidate rectangle (perhaps favouring rectangles with their centre 
near the centre of the X-ray image). 

- Orientation: The orientation of the candidate rectangle within the image (perhaps favouring rectangles parallel with 
the edge of the X-ray image). 

10 

Region statistics: 

- Exposure: The mean, maximum, minimum (and the like) statistical measures derived from the (candidate) signal 
and shadow regions of the X-ray image (perhaps favouring tow-exposure grey-levels in the shadow, versus signal 
regions). 

Homogeneity: The variance (and the like) statistics derived from the (candidate) signal and shadow regions of the 
X-ray image (perhaps favouring low-variance grey-levels in the shadow, versus signal regions). 

Topology: 

20 

Visibility: The visibility of the candidate signal/shadow boundary, i.e. the degree to which the candidate signal/ 
shadow boundary is wholly contained within the image (perhaps favouring complete, rather than clipped boundaries 
- eg as in figure 1d rather than fig. 1e). 
Completeness: The number of visible rectangle sides. 

25 

Hypotheses consisting of just two 'sides' (which are approximately orthogonal or parallel), can be considered as partially 
visible rectangles. The scoring may, for instance, favour hypotheses with more, over fewer, sides. 

The scalar result of each test may provide for the rejection of a hypothesis by assuming an appropriately small 
positive or large negative value, effectively denying that hypothesis any chance of highest ranking. 
30 Step v) is just the selection of that hypothesis which is ranked the highest of those surviving the high-level tests 

of step iv). 

The method of the present invention can be applied for recognising an irradiation field in cases where at the 
exposure step an irradiation field stop has been used. 

It will be clear that the invention is not limited to this particular application. It can for example also be used when 
35 the borders of irradiation fields have to be determined in an application wherein at least two images are recorded in a 
juxtaposed relation, a mode often referred to as 'split screen'. 

Brief description of drawings 

40 The-method of the present invention as well as particular embodiments hereof will be explained hereinafter with 
reference to the accompanying drawings wherein 

Figure 1 shows examples of location, orientation, and visibility of signal/shadow boundary, 

45 Figure 1a shows an image with no collimation shadow regions, 

Figure 1b shows an image wherein two sides of the rectangle are visible, (1) is the signal region, (2) are the 
shadow regions, 

Figure 1c shows an image with 3 sides of the rectangle visible, 
Figure 1d shows an image with all 4 sides of the rectangle visible, 
50 Figure 1 e shows an image wherein 4 sides are partially visible, 

Figure 1f shows an image with 1 side visible, 

Figure 2 is an example of Canny Edge Detection, 

Figure 3 is an example for line clustering. Lines l b and l c may be clustered. Similarly, lines l e , If, and l g may be 
55 clustered. 

Figure 4 is an example of edge thinning, 

Figure 5 are examples of vertex classifications for the central pixel in each 3x3 array, 
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Fig.5a:isolated point 
Fig.5b:terminal, 
Fig.5c:edge-point, 
Fig.5d:T-junctlon, 
5 Fig.5e:X-junction. 

Figure 6 is an example of linear extrapolation. Points P a and P b are the limits of line \ B . 
Figure 7 are examples of closest approach of two lines, 
Figure 8 is an example showing overlapping, and non-overlapping lines. 
10 Figure 9 is an example showing edge image, vertex image edge-points, and projection of line-clusters, 

Figure 10 is a possible hypotheses for the example of figure 9, using clusters L^, L^, and L c , 
Figure 11 is a possible search tree for example of figures 9 and 10, the letters denote the labelling in figure 10, 
Figure 12 is an exposure and readout apparatus. 

*5 The following is a detailed description of each of the individual method steps of the present invention. 

Step i) Extraction of low-level primitives 

The preferred low-level primitives are lines. They are derived from the X-ray image in two stages: edge-detection 
20 and segmentation. 

1.1 Edge Detection 

The purpose of this stage is to produce an edge-image: see for example figure 9a. Our system incorporates a 
25 Canny edge-detection process. 

The following are definitions of some terms used in connection with a Canny edge detection process: 
Canny edge-detection: A particular process for edge-detection comprising optional smoothing, followed by pro- 
duction of an edge-image wherein a pixel is labelled as an edge-pixel if and only if the Sobel squared magnitude at 
that point is a local maximum along the (one-dimensional) trajectory which is normal to the Sobel orientation. 
30 A Canny edge-image {E(r)J is formed by non-maximal suppression of the Sobel squared magnitude in the direction 

of the orientation of the image gradient, i.e. normal to the Sobel orientation: 

E(£) = 1 if G(r) > G(r+R q ) > G(r-R q ) 

35 

= 0 otherwise 

where 

40 R =(1,0)if IIQ(r)l - ji/2 ISn/8 

= (1,-1) if IIQ(r)+7i/4l- jc/2 I < 7t/B 
= (1,1) if IQ(r)+7i/4l<*/8 
= (0,1) if lQ(r)l<7t/8 

Let a mean a is bold and vectors r, R q represent the two-dimensional coordinates on the image plane (r= (i,j) say). 
so For an illustration, see figure 2. 

Smoothed image: An intrinsic image derived from the X-ray image wherein the high spatial frequencies are sup- 
pressed relative to the low spatial frequencies. E.g. the result of convolution of the X-ray image with a (two-dimensional) 
Gaussian. 

Sobel squared magnitude: The pixel-value G(i,j) representing the square of the gradient-magnitude of the image 
55 (X(i,j)} as determined by the following formula: 

G(ij) = (X(i+1 ,j) - X(i-1 ,j)) 2 + (X(i, j + 1 ) . X(i,j-1 )) 2 
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Sobel orientation: The pixel-value Q(i,j) representing the orientation of the tangent to the iso-gradient-magnitude 
of the image {X(i,j)}, as determined by the following formula: 

5 Q(i,j) = atan2(X(i+1 ,j) - X(i-1 X(i.j+1 ) - X(i,j-1 )) 

where atan2(y,x) is the arctangent C language function defined in (-71,71). 

The edge-detection process incorporated in our system comprises four stages: 

10 1.1.1 Subsampling of the X-ray image 

The subsampling reduces the number of pixels of an image. For example: an original image obtained by reading 
an irradiation image stored in a photostimulable phosphor screen (cfr. figure 12) comprising approximately 5 million 
pixels, is reduced to an image of about 20000 12-bit pixels, preserving the aspect ratio of the original image. 

15 

1.1.2 Canny Edge-Detection 

The subsampled image is passed to a Canny edge-detector to produce images: Sobel-squared magnitude {G(i, 
j)}, Sobel-orientation {Q(i j)}, and edge-image {E(i,j)}. The Gaussian smoothing is via a symmetric 5x5 convolution 
20 mask with weights calculated to effect a standard deviation of 0.35 pixels. 

1 .1 .3 Elimination of Weak Edges 

Edge-points with a low gradient magnitude indicate a weak edge. Studies have shown that the signal/shadow 
25 boundary is generally of high contrast, giving a large gradient magnitude. Very weak edges are therefore unlikely to 
form part of the signal/shadow boundary. Those with a Sobel-squared-magnitude below a threshold of 2000 are re- 
moved from {E(i,j)}, and correspondingly from {Q(i,j)}, and {G(i,j)J. 

1.1.4 Angle of Locus 

30 

The orientation {O(i.j)}, of the locus of edge-points in the edge-image {E(ij)}, is calculated in three stages: edge- 
thinning, connectivity analysis, linear regression. 

1.1.4.1 Thinning 

35 

The edge-image is destructively thinned. The motivation is threefold: 

1 . The parameters of a regressed straight line, including the goodness-of-fit, are more accurately determined from 
the spatial distribution of edge-points. 
40 2. The various types of junction (TERMINAL, T-JUNCTION, ...) are more easily detected. This is required for corner- 

detection (below). 

3. There then exists a close correspondence between the length of a straight line, and the number of (connected) 
edge-pixels of which it is made up. This is required for line-statistics (above). 

45 Thinning is achieved by removing edge-pixels which do not affect the local (3x3) connectivity of other edge-pixels. 

See figure 4 for an example. 

1 .1 .4.2 Connectivity Analysis 

so Connected edges are uniquely labelled by recursive edge-walking to create a label image {L(i,j)}. All edge-pixels 

mutually connected share a common, but otherwise unique, label. 

1.1.4.3 Linear regression 

55 in this connection the following definitions are first given: 

Linear regression: The process by which an edge can be optimally parameterised according to the spatial location 
of its constituent edge-pixels. E.g. the determination of the parameters of a line which are the radius to the normal 
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IS 



20 



(from the rigin), and its orientation r lative to the X-axis, by a least squares (of perp ndicular distances) fit. 

W ighted linear regression: The proc ss by which an edge can b optimally parametehsed according to both the 
spatial location of its constituent edg -pixels, and th degree of certainty as to th ir correct classification (labelling). 
E.g. the determination of the parameters of a line which are the radius F^r to the normal (from th origin), and 
its orientation O^r relative to the x-axis, by the minimisation of a weighted least squares (of perpendicular dis- 
tances) fit, the weighting being the Sobel squared magnitude. 

Let p k be the positions of the edge-pixels (indexed by k) which comprise the line: 
and let the positional covariance matrix be P. 



F» r (x,x)^ r <x,y) 

I r <y,x)^ p(y/y) J 



given by T^^7^^^N k - 2^(^2^( b >w^ for a,b in {x,y} then Qwlr = 1/2 atan2(-2r< x 'y>, r^-rM) and 
= (sin(Q WLR )Z k p k ( x )w k + costQwLRjZ^MwkJ/ZkWk The eigenvalues \ + \ of r are the sum of square errors (i.e. pro- 
25 portional to the variances) in each of the principal directions ie along, and normal to, the best fit straight line: 

2X ± =Tr(r)± (Tr(D 2 -4Det(r)) 1/2 

30 in the described embodiment of the present invention the orientation of the locus at each point is calculated by 

linear regression of connected edge-points in a local neighbourhood (nominally 7x7). The connectivity is guaranteed 
by including only those pixels which, in {L(i,j)), share a common label. The positional-covariance matrix gives the 
orientation of the maximum variance of the data which is also the orientation of the least-squares fit straight line with 
two degrees of freedom. 

35 With reference to the definition of weighted linear regression, with w^1 : 0(i,j)=Q WLR is the orientation of the locus 

of edge-points without regard to the sense of the original edge. The ambiguity in the sense of the locus is resolved by 
borrowing from the Sobel orientation Q(r) : the direction of the tangent to the Sobel squared magnitude iso-contour is 
always close to the direction of the tangent to the locus of edge-points as determined by the Canny edge-detector. 
Therefore the orientation is adjusted to minimise the angular separation between 0(i,j) and Q(i,j): 

40 

Min {0(i,j)-Q(iJ) + kn}^k min 
k in {-3..3J 

45 

O(iJ) :=0(i,j)-Q(i,j) + k min 7c+2njc 



where n is chosen so that the result is in [0,2n]. 

so 

1 .2 Segmentation 

Line primitives are built from the remaining edges. The chains of edge-pixels are broken into piece-wise straight 
sections (figure 9b), and the pixels from each section assembled into a line structure. 
55 There are four steps to this procedure: 
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1.2.1 Vertex Analysis 

An intrinsic image {V(i,j)) is constructed with labels for the classification of edge-pixels in {E(i,j)}. According to th 
3x3 connectivity of edge-pix Is in {E(i,j)}, one of the following labels are awarded to the corresponding point in {V(i,j)}: 
5 ISOLATED-POINT, TERMINAL, EDGE-POINT, T-JUNCTION, X-JUNCTION. An example of each is given in figure 5. 

1 .2.2 Comer Detection 

Comers are detected by a refinement of the EDGE-POINT labels in {V(i,j)|. There are two stages: 

w 

1.2.2.1 Curve detection 

An intrinsic image (U(i,j)} is calculated with values assigned to each pixel which corresponds to an EDGE-POINT 
labelling in {V(i,j)}. The pixel values in {U(i,j)} measure the local conformity of edge-pixels in {E(i,j)}, in a 7x7 neigh- 
bourhood, to a straight-line. The positional covariance matrix is calculated as described in 1 .1 .4.3 with all w^1 . 

The conformity, (to a straight line) is taken to be the ratio of major to minor eigenvalues of the positional covariance 
matrix, ie Pixels previously labelled as EDGE-POINTs in {V(i,j)} are then relabelled as CURVES, if the conformity 
is less than a threshold value: 

If V(i,j)=EDGE-POINT , and KK > 0-04, then V(i,j) = CURVE. 

20 

1.2.2.2 Non-maximal supression 

The pixels in {U(i,j)} are identified as comers if they are both CURVE points in {V(i,j)} (i.e. having a value above 
some threshold), and they are locally maximum in a 3x3 neighbourhood. The corresponding pixels in {V(i,j)} are rela- 
ys belled as CORNERS. 

1 .2.2 Extraction of connected edges 

Those pixels in {V(i,j)} still retaining an EDGE-POINT classification are candidates for incorporation into the low- 
30 level primitives. Such pixels are extracted from |V(i,j)}, into {S(i,j)} say, which is the linear segmentation as an intrinsic 
image: each connected chain of pixels is given a unique label in {S(i,j)}. Each connected set of pixels is then gathered 
into a line structure which is the low-level primitive. 

1.2.2.1 Rejection of short lines 

35 

Studies have shown that sides of the signal/shadow boundary generally have at least some portion of their length 
for which the edge-detector furnishes an unbroken chain of edge-pixels. In order to avoid false positives therefore, an 
implicit rule is that a side cannot be composed entirely of very short edge-chains, but must have at least one line (in 
the line-group) which is of substantial length. Prior to segmentation and extraction of the line primitives therefore, chains 
40 of edge-pixels shorter than some threshold length are rejected (from (S(i,j)}) at this point. By this means, the number 
of low-level primitives is reduced without compromising the accuracy: short lines may form part of the final hypothesis 
following data-retrieval. 

1 .2.3 Weighted linear regression 

46 

A weighted linear regression is performed on each of the lines, by which means the optimal parameters (normal 
radius to origin, and orientation), and conformity are recorded. The weights are the Sobel-squared magnitudes of the 
contributing edge-pixels, as a fraction of the maximum Sobel-squared magnitude: 

Mi w k (i,j) = G(i,jKMax{G(u,v)} V i,j 

Vu,v 

55 1 .2.3.1 Rejection of poor lines 

Lines are no longer 'local' entities. The test on conformity to a straight line is now repeated with a less tolerant 
threshold in accordance with the extra precision afforded by the (arbitrarily large) number of edge-pixels in the line: 
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1.2.4 Linear xtrapolaiion 

The surviving Jines are extrapolated from each end, until they intersect either the edge of the imag , or an edge- 
point which does not belong to that line (as determined from the labels in {S (i.j)}). The 'limits' so determined, are saved 
5 in the line structure. A distinction is made between each of the two limits which correspond to the two ends of the line. 
This distinction respects the sense of the line: see figure 6. 

1 .2.5 Data-retrieval 

10 Edge-points are now recovered from {E(i, j)} which were not classified as EDGE-POINTs in {V(i,j)} (perhaps because 

they were reclassified as CORNERS, CURVEs or JUNCTIONS etc). For each line, two rectangular areas are scanned 
in {E(i,j)J each of which is centred on (aligned with) the projection of the line as given by the regressed parameters 
r wlr and Qwlr T" 9 rectangles extend from the line end-points, to the limits as determined by extrapolation. The 
width of the rectangle is 5 pixels. All edge-points in {E(i,j)} which were not already attributed to the line are now included. 

15 

Step ii) Extraction of intermediate-level primitives 
2.1 Linear agglomerative clustering 

20 In this connection the following definition are given: 

"Linear agglomerative clustering" generally is an iterative process for grouping together lines which are deemed 
(by some criteria) to be sufficiently similar. E.g. a pair- wise agglomerative clustering process based on an analysis of 
variance test, and broadly as described in "Unsupervised Estimation of Polynomial Approximations To Smooth Surfaces 
in Images or Range Data" by J.F. Silverman and D.B. Cooper, Proc. IEEE Cont. on Automation and Robotics, 19B6, 

2S page 299-304. and defined hereinbelow, but with the mergable entities being lines rather than regions, and the con- 
straint that two mergable regions be touching replaced by a constraint that two lines (or line-clusters) be linearly co- 
herent. 

"Analysis of variance" is a statistical method for comparing (the similarity of) two sample distributions under the 
assumption that the underlying distribution is Gaussian; but the variance is unknown. The method determines an F- 
30 statistic, which is the ratio of the two sample variances, correctly normalised to take account of the true degrees of 
freedom. 

When the test is to determine the likelihood of two distributions having the same parent distribution, then the ratio 
is between the unexplained and the explained variances. When applied to pair-wise clustering, the unexplained vari- 
ance is the extra that is incurred by the proposed merge: 

35 

<°u 2> = (* m 2 - *a 2 * * b 2 )/lv m " v a -v b l 

where a,b and m denote the separate, and merged entities respectively; x and v are the chi-squared variates and 
40 corresponding degrees of freedom. The explained variance is that prior to the proposed merge: 

«* e 2 > = (X a 2 + X b 2 V( v a + v b) 

45 and 

F(v u ,v e ) = <o u 2 >/<a e 2 > 

so is the F statistic, having degrees of freedom 

v = Iv - V - Vul 

55 and 

v e = v a + v b 
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The F statistic is used to index the cumulative distribution function Q(Frv a ,v b ) which gives the probability of that 
value or greater occuring by chance. If the un xplained variance (as a fraction of the xplained) is low, then F is small, 
and Q(F!v a ,v b l) is larg . In that case, the entities may be safely clustered. 

In the described mbodiment of the pres nt invention lin s underg linear agglomerative clustering according to 
5 a test of their linear coherency. For all line or line-cluster pairs, if the pair passes the test of linear coherency, an F 
statistic based on an analysis of variance, is computed. The pair with the highest probability Q(Frv a ,v b ) is merged. The 
parameters and limits of the merged line-cluster are recomputed using weighted linear regression. The whole cycle is 
repeated until no line-clusters pass the linear coherency test. An example of a set of line-clusters, following linear 
aggbmerative clustering, is given in figure 9c. 

10 

2.1.1 Linear Coherency 
Definitions: 

is Linearly coherent: Two lines or line-clusters are said to be linearly coherent if they could potentially form part of 

the same side of the rectangular signal/shadow boundary. 

Linear agglomerative clustering: an iterative process for grouping together lines which are deemed (by some criteria) 
to be sufficiently similar. E.g. a pair-wise agglomerative clustering process based on an analysis of variance test, and 
broadly as described in (Silverman & Cooper), but with the mergable entities being lines rather than regions, and the 
20 constraint that two mergable regions be touching replaced by a constraint that two lines (or line-clusters) be linearly 
coherent. 

In the described embodiment of the present invention the linear coherency test is applied to each line or line-cluster 
pair during linear agglomerative clustering. For two line or line-clusters to be viable candidates for clustering: 

25 1 . They must be approximately coiinear. 

2. Their extrapolations must overlap or touch. 

2.1.1.1 Colinearity 

30 Two line-clusters are determined as coiinear if: 

1 . The line-clusters are approximately parallel: 

35 Min{Q (1) WLR -Q (2) wtR .2en}<15» 

ein {-1,0,1} 
and 

2. The distance of closest approach between the two line-clusters is less than some threshold: 
40 d c < 4 pixels where the d c are the distances depicted in figure 7. 

2.1.1.2 Overlapping lines 

Two line-clusters are said to have overlapping extrapolations if either has a perpendicular projection onto the other: 
45 see figure 8. 

2.1.2 Analysis Of Variance 

The F-statistic based upon the analysis of variance is applied to test the likelihood that two sample distributions 
50 come from the same parent distribution. With reference to the definition of analysis of variance, and weighted linear 
regression, the chi-squared variates are the sum of square errors of the best fit straight line: 




where the line-cluster index j is any of a.b.m. The corresponding degrees of freedom are: 
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where N a and N b are the numbers of edge-pixels In the lines (of labels a and b respectively). 
10 Step hi) Building hypotheses 

3.1. Search Tree 

A search tree is constructed of all legal combinations of line-clusters : see figure 10a-h. The tree starts with a root 
is node which corresponds to the null hypothesis that there is no collimation material present, and therefore that the 
signal/shadow boundary is absent. The addition of a new branch & node to the tree represents the addition of a further 
line-cluster to the hypothesis inherited from the parent node : each new node is associated with a new line-cluster. 
Therefore the return path from a particular node to the root determines the line-clusters in the hypothesis associated 
with that node. 

20 The search tree has at most five levels, corresponding to the null hypothesis, and four levels for each of the four 

sides of the rectangular signal/shadow boundary. Each node represents a hypothesis in its own right. 

As each node is added, the intermediate-level tests are applied to test the legality of the hypothesis. If the results 
are negative, then the node is flagged as illegal, excluded from further consideration, and the creation of offspring 
prohibited. Assuming all the hypotheses shown are legal, the example of figure 10 would result in a search tree in its 

25 final state as given in figure 11 . 

3.2. Intermediate-level tests 

The following must be satisfied by each hypothesis during the construction of the search tree : 

30 

3.2.1. Orthogonality 



35 



40 



45 



All sides in a hypothesis must be approximately mutually orthogonal. Following addition of a new line-cluster at 
node n, this test is applied as : 

Min {l(Q (n> wlr-Q^ WLR -kn/2l}< 0.13 
kin {-4. -1,1 ...4} 

where the Q are in radians, and where p ranges over ail nodes but n in the return path. If the new node fails this test, 
it is declared illegal. 

3.2.2. Closure 

All sides in a hypothesis must meet, at each end, either another side, or the edge of the image. Following the 
addition of a new line-cluster at node n, this test is applied by examining the limits of the new line-cluster at node n, 
and other line-clusters (at nodes p) in the return path which form an adjacent side of the rectangle. 



so Step iv) Testing and costing hypotheses 

Once the hypothesis has been built, the intersection points for each of the sides of the candidate signal/shadow 
boundary are computed. These intersections may be either with the edge of the image, or with another side; see figure 
9. The high level rules must then be satisfied in order that the hypothesis is verified. 

55 

4. 1 . Definition of terms and symbols 

Image dimensions - Let th number of pixels in the image array be Imagewidth and Imageheight in the x and y 
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directions respectively. 

Numb r of sides - L t the number of (visible) sides in the hypothesis be N (O £ N £ 4). 

Intersection points - Let the intersection points of the kth side be (cc^) and (^,6 k ) for the int rsections nearest 
the lower, and upper limits, l k and u k respectively. 
Interval - Let the interval of the kth side b 



RECTANGLE NOT SMALL - Let the signal/shadowboundary not be small if A k > = 10 for all k in {1..N} 

Accumulated evidence - Let the accumulated evidence for the kth side be W k = ZjW k j where the sum is over all 
edge-points contributing to the kth side. Let the accumulated evidence for the kth interval be W k = XjW ki where the 
sum is over all edge-points contributing in the interval of the kth side. 

EVIDENCE WELL USED - Let the amount of 'unused' evidence be small if: WyW k > = 0.4 for all k in {1 ..N} 

Exposure statistics - Let max k and mean k be the maximum and mean exposures (in the X-ray image) in the half- 
plane exterior to the kth side of the rectangle. Let max^ and mean^ be the maximum and mean exposures in the 
interior (signal) part of the rectangle. 

STRONG CONTRAST BOUNDARY - Let the boundary have strong contrast if : (maXk-MinfXti^jyfmaXi^lvlinlXfi, 
j)}) < 0.65 and mean^mean^ < 1 .05 for all k in {1 ..N} 

VERY STRONG CONTRAST BOUNDARY - Let the boundary have very strong contrast if : 
(max k -Min{X(i,j)})/(max jn -Min{X(i,j)}) < 0.3 and 
meani/meanjn < 1 .05 for all k in {1 ..N} 

Line support - let pix k be the number of pixels in the interval of the kth side, for which some evidence exists, ie for 
which w^O. 

PERSISTENT BOUNDARY - Let the boundary be persistent if : 
(pix k + 2)/A k > 0.8 for all k in {1 ..N} 

VERY PERSISTENT BOUNDARY - Let the boundary be very persistent if : (pix k + 2)/A k > 0.88 for all k in {1 ..N} 
BOUNDARY HAS TWO ALIGNED SYMMETRIC SIDES - Let the boundary have two aligned symmetric sides if : 



A k = Max { IcyYfcUVM + 1 - 



N = 2 and 




f 



and 




RjJ[ R l/lmagewidth<0.1 



or 



N = 2 
and 




and 




FT 



{ (2) 



I / Imagewidth < 0.2 
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where th labels K = 1 ,2,3,4 corr spond to those lines having angl s with respect to the x axis, closest to id2 t n, 3nJ 
2, and 2n radians. 

BOUNDARY HAS STRONG GESTALT - The boundary has strong gestalt if : 
N = 4 or BOUNDARY HAS TWO ALIGNED SYMMETRIC SIDES 

BOUNDARY HAS FOUR ALIGNED AND ROUGHLY SYMMETRIC SIDES - Tru if : 
N = 4 
and 



10 



Min f 

I 

n in {-1,0,1} [ 



(K) , 
Q WLR - Krc/2 - 2n* 1/ < 0.087 



VK in {1,4} 



15 



20 



25 



30 



and 



and 



I RjJ^ - R^ R I / Imagewidth < 0.2 



I R|^ R - R^ R I / Imagewidth < 0.2 



BOUNDARY HAS TWO LONG ASPECT SYMMETRIC SIDES - True if : 
N = 2 

and Imageheight £ Imagewidth 
and 



35 



Min 

n in {-1,0,1} 



I Q (K> - KK/2 - 2nn I 



< 0.087 / VK in {1,3} 



and 



40 



I R^ LR - R^ R I / Imagewidth < 0.1 



orN = 2 

and Imagewidth < Imageheight 
45 and 



50 



Min [ 

. i 

n in {-1,0,1} I 



I Q (K> - K7C/2 - 2n* I 



< 0.087 



VK in {2,4} 



and 



55 



I RjJ[ R - R^ R I / Imagewidth < 0.1 



BOUNDARY HAS VERY STRONG GESTALT - True if : 
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BOUNDARY HAS FOUR ALIGNED AND ROUGHLY SYMMETRIC SIDES or BOUNDARY HAS TWO LONG ASPECT 
SYMMETRIC SIDES 

4.2. Verification rules 

5 

A hypothesis is accepted if the following composite rule is satisfied : 
RECTANGLE NOT SMALL and EVIDENCE WELL USED and (BOUNDARY HAS VERY STRONG CONTRAST and 
BOUNDARY IS PERSISTENT or BOUNDARY HAS STRONG CONTRAST and BOUNDARY IS VERY PERSISTENT 
or BOUNDARY HAS VERY STRONG GESTALT and BOUNDARY IS PERSISTENT or BOUNDARY HAS VERY 
io STRONG GESTALT and BOUNDARY HAS STRONG CONTRAST or BOUNDARY HAS STRONG GESTALT and 
BOUNDARY IS VERY PERSISTENT) 

4.3. Costing hypotheses 

*5 |f a hypothesis is rejected, it is awarded (effectively) an infinitely large cost. If the hypothesis is accepted however, 

it is awarded a cost which is determined as described below : 

Let p ou! = 0.1 be the fraction of pixels on the collimation/signal boundary which do not have a corresponding edge- 
pixel as determined by the Canny edge-detector. 

Let p in = 0.1 be the fraction of pixels which are co-linear with, but not on, a side of the collimation/signal boundary, 
20 yet which are labelled as edge-pixels by the Canny edge-detector. 

Then the probability P^ of having k or more edge-pixels missing on a collimation/signal boundary of total interval 
(between intersects) n pixels is : 



25 




p s q n-s . Q(r|u lr n 2 ) 



30 

where q = 1-P in ; w 1 =2(n-k+1); x> 2 =2k; F^kO-p^fn-k+l). 

The same formula applies to the probability P ou! , of there being k or more edge-pixels which are colinear with, but 
not on, sides of the collimation/signal boundary, of total interval between intersects n, and total interval between limits 
of L, provided : 

35 

* = 1 -Pour »i = 2(L-n-k + 1); u 2 = 2k; F = k(1-p out y(L-n-k+1). 
The cost awarded is then : 

40 

4 - N + Min{0.999 , -log(P ift P out )/100}. 

This function constrains the contribution from missing, and extra, edge-pixels to lie in [0,0.999]. The main part of 
45 the cost is determined by the number of sides : 0 for 4 sides, 1 for 3 sides, etc. The implicit rules area that : 

1 . Validated hypotheses are unconditionally cheaper for having more sides. 

2. Given two hypothesis of equal number of sides, the cheapest is that which best explains the data. 

50 Step v) Choosing the best hypothesis 

The solution search tree is maintained by a branch and bound function. All possible hypotheses (potentially legal 
combinations of line-clusters) are considered. The hypothesis chosen by the system is that with the lowest cost, as 
awarded at step iv). If no other hypothesis is validated, the null hypothesis that there is not signal/shadow boundary 
55 (and therefore no shadow region) is always valid (with a cost of 4 units). 

Figure 12 shows an embodiment of an exposure and read-out apparatus. A radiation image of an object was 
recorded on a photostimulable phosphor screen (6) by exposing (3) said screen to X-rays transmitted through the 
object. The irradiation field emitted by the X-ray source was limited by using a opaque X-ray shield (4). 
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The radiation image stored in the photostimulable phosphor scr en was read-out by scanning the phosphor screen 
with stimulating rays emitted by a laser (9). The stimulating rays w re deflected according to the main scanning dir ction 
by means of galvanometric deflection d vice (8). Th secondary scanning motion was p rform d by transporting the 
phosphor scr en in a direction indicated by arrow (7) perpendicular to the main scanning direction. The stimulated 
5 emission was directed by means of light collector (1 0) onto a photomultiplier (11 ) where it was conv rted into an electric 
signal representation. The analog output signal is then digitised (not shown). 

This signal can then be subjected to a processing method according to the present invention for the purpose of deter- 
mining the location of the signal -shadow boundary and hence also to the determination of the useful image region. 
Contrast enhancing processing can then be limited to the diagnostically useful image part. 

10 

Claims 

1. A method for the determination of the location of the signal/shadow boundary in an X-ray image represented in a 
is digital signal representation comprising the steps of : 

i) Extracting low-level primitives from the X-ray image {X(i,j)}» 

ii) Forming a reduced number of intermediate-level primitives from the low-level primitives, 

iii) Building hypotheses being possible solutions regarding the location of the signal-shadow boundary, from 
20 combinations of intermediate-level primitives, during which each combination is subject to intermediate-level 

tests, whereupon partial or complete hypotheses are rejected or accepted, 

iv) Performing high-level verification tests on each hypothesis, whereupon hypotheses are rejected, or accept- 
ed at some cost, said cost defining a ranking of the acceptability of a hypothesis, 

v) Selecting the hypothesis with the least cost. 

25 

2. A method according to claim 1 wherein said low level primitives are lines. 

3. A method according to claim 2 wherein said lines are extracted by performing on the digital signal representation 
the following steps: 

30 

subsampling, 

performing edge detection on said subsampled image to produce the following intrinsic images: 

squared local gradient magnitude {G(i.j)}, 

35 - local edge orientation {O(i.j)}, 

edge image {E(i,j)}, in which all pixels having a squared gradient magnitude below a predetermined thresh- 
old value are set to zero, and in which all edge pixels which do not affect the local connectivity of other 
edge pixels have been removed by a thinning operation, 
- label image {L(i,j)J, in which all mutually connected edge pixels from the thinned edge image share a 

40 common, but otherwise unique label, 

local edge orientation image, obtained by linear regression of connected edge pixels in a local neighbour- 
hood, 

vertex image {V(i,j)}, in which all edge pixels are identified to be one of the following vertex types: isolated 
point, terminal, edge point, T-junction, X-junction, curve point, or comer, 
45 - segment image {S(i,j)}, in which all edge pixels of vertex type 'edge point 1 , as obtained from {V(i,j)}, are 

assigned the label of the corresponding pixel in (L(i,j)}, and the other pixels are zero, 

extracting from {S(i,j)} pixels with identical labels, 

creating line data structures for the purpose of representing lines, said data structures comprising fields for 
so representing other line-specific data such as required by subsequent method steps and tests, 

storing each set of coordinates of identically labelled pixels extracted from S{(i,j)} in a line data structure, 
estimating for each of the line data structures the length of the corresponding line, 

eliminating line data structures corresponding to all lines having a length smaller than some predetermined 
value from said line data structures, 
55 - determining the normal radius and orientation of each line by linear regression and recording those values in 

said line data structures as attributes of the corresponding line, 

eliminating all line data structures for which the conformity to a straight line as determined by linear regression 
is below some predetermined value, 
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lin arly extrapolating each non-eliminated lin from each nd, until it intersects either the image border or with 
an edg point as determined by {S(i,j)J which does not belong to that line, and recording in th line data 
structures the positions of both int rsection points as attributes of the corresponding lin . 

5 4. A method according to claim 3 wher in said identification of the vertex type CORNER comprises 

- curvature determination of connected chains of edge-pixels excluding pixels which are junctions or terminators, 
local non-maximum suppression of the curvature values, 

- threshold of the remaining non-zero curvature values, accepting only those points as corners whose value is 
io above some threshold. 

5. A method according to claim 3 wherein said identification of the vertex type CURVE comprises 

curvature determination of connected chains of edge-pixels excluding pixels which are junctions or terminators, 
is - threshold of the curvature values, accepting only those points as curves whose value is above some threshold. 

6. A method according to claim 3 wherein weighted linear regression comprises: 

- calculation of weighted positional moments, upto second order, of edge-pixels in each line data structure, 
20 wherein the weighting of each pixel contribution at (i,j) is a functbn of the gradient magnitude G(i,j), 

use of the weighted moments to determine the parameters of the best fit straight line so as to minimise the 
weighted sum of squared errors (of pixel positions to the parameterised line). 

7. A method according to claim 2 wherein said intermediate level primitives are line groups obtained by a linear 
25 agglomerative clustering process. . 

8. A method according to claim 7 wherein in the linear agglomerative clustering process any two lines or two line 
groups, or a line and a line group are merged into a single line group if: 

30 - their difference of orientation is less than some predetermined value, 

their distance of closest approach is less than some predetermined value, 
they have overlapping extrapolations 
• the outcome of the F-statistic based upon the analysis of variance is such that the two considered pixel dis- 
tributions most likely belong to a single linear pixel distribution which constitutes the merged cluster. 

35 

9. A method according to claim 7, where candidate hypotheses for the signal/shadow boundary are formed from 
every possible line-group, pair of line-groups, triplet of line groups, and quartet of line-groups, including, in addition, 
the hypothesis that there is no signal/shadow boundary. 

40 10. A method according to claim 7, wherein said high-level verification tests for the acceptability of a hypothetical 
signal/shadow boundary, comprise any logical combinations of the following tests for 

boundary closure, for which test to be passed, for the hypothesis under consideration, the extrapolated line- 
groups must each terminate on, or intersect with, another line-group in the same hypothesis, or the edge of 
45 the image, 

boundary shape, for which test to be passed, for the hypothesis under consideration, is constrained to be 
rectangular, or nearly rectangular, or in the case of there being less than four sides in the hypothesis, the 
visible sides are constrained to be consistent with a partially visible rectangular or nearly rectangular shape, 
signal/shadow contrast, for which test to be passed, for the hypothesis under consideration, is constrained to 

so be largely consistent with there being a lower level of exposure in the region of {X(i,j)} corresponding to the 

shadow region, than the level of exposure in the region of {X(i,j)} corresponding to the signal region, 
boundary location, for which test to be passed, for the hypothesis under consideration, the centre of the image 
must be either contained within the signal region, or must be close to a signal/shadow boundary, 
boundary alignment, for which test to be passed, for the hypothesis under consideration, each visible side of 

55 th rectangle must be parallel to a side of the image, 

- boundary symmetry, for which test to be passed, for the hypothesis under consid ration, there must be two 
or four visible sides, which are aligned (as above), and approximately equidistant from the edge of the image, 
boundary completion, for which test to be passed, for the hypothesis under consideration, the fraction of the 
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boundary perimeter length, for which ther exists evidence from the low-level primitlv s, must be gr ater than 
some threshold. 

1 1 . A method according to claim 2 wherein a cost is assigned to each hypothesis, which cost is low st if the hypothesis 
s is that of a rectangle with four sides visible, and which cost increas s as the number of visible sides of a rectangle 

in a hypothesis decreases, and the cost is highest if the hypothesis has no visible sides. 

12. A method according to claim 11, for further refining the cost of each hypothesis, wherein the cost is increased 
according to the fraction of the signal/shadow boundary which is not supported by evidence from the low-level 

10 primitives. 

13. A method according to any of the preceding claims wherein image processing performed on said digital signal 
representation of said X-ray image is limited to pixel values of said image within the signal/shadow boundary. 

15 14. In a method of reading a radiation image stored in a photostimulable phosphor screen by scanning said screen 
with stimulating irradiation, detecting light emitted upon stimulation and converting the detected light into a digital 
signal representation a method of determining a signal/shadow boundary according to any of the preceding claims. 



20 Patentanspruche 

1. Ein Verfahren zur Bestimmung des Ortes der Signal-Schatten-Grenze in einem in einer digitalen Signaldarstellung 
dargestellten Rontgenbild, wobei das Verfahren folgende Schritte umfaSt: 

25 j) Ableiten von Grundelementen niedriger Stufe aus dem Rontgenbild {X(i,j)}, 

ii) Bilden einer verringerten Anzahl von Grundelementen der Zwischenstufe aus den Grundelementen niedri- 
ger Stufe, 

iii) Bilden von Hypothesen, die mogliche Losungen hinsichtlich des Ortes der Signal-Schatten-Grenze sind, 
aus Kombinationen von Grundelementen der Zwischenstufe, wahrend dessen jede Kombination Zwischen- 

30 stufentests unterworf en wird, worauf hin teilweise Oder vollstandige Hypothesen verworf en oder akzeptiert wer- 

den, 

iv) jede Hypothese hochstufigen Uberprufungstests unterziehen, woraufhin Hypothesen verworfen Oder unter 
gewissen Kosten akzeptiert werden, wobei die Kosten eine rangmaBige Einstuf ung der Akzeptierbarkeit einer 
Hypothese definieren, 

35 v) Auswahl der Hypothese mit den geringsten Kosten. 

2. Verfahren nach Anspruch 1 , wobei es sich bei den Grundelementen niedriger Stufe urn Linien handelt. 

3. Verfahren nach Anpruch 2, bei dem die Linien abgeleitet werden, indem an der digitalen Signaldarstellung die 
40 folgenden Schritte durchgefuhrt werden: 

Unterabtastung, 

Durchf uhrung einer Randerfassung an dem unterabgetasteten Bild, urn die folgenden intrinsischen Bilder her- 
zustellen: 

45 

~ quadierte GroBe {G(i,j)} des ortlichen Gradienten, 

- ortliche Randausrichtung {O(i.j)}, 

- Randbild {E(i,j)}, in dem alle Bildpunkte mit einer quadrierten Gradientengrd&e unterhalb eines vorbe- 
stimmten Schwellwerts auf null gesetzt werden und bei dem alle Randbildpunkte, die sich nicht auf die 

so ortliche Verknupfbarkeit anderer Randbildpunkte auswirken, durch einen AusdOnnungsvorgang entfemt 

worden sind, 

Markenbild {L(i, j)}, bei dem alle miteinander verbundenen Randbildpunkte von dem ausgedunnten Rand- 
bild eine gemeinsame, aber ansonsten eindeutige Marke gemeinsam haben, 

ortliches Randausrichtungsbild, das durch lineare Regression verbundener Randbildpunkte in einer ort- 
ss lichen Umgebung gewonnen worden ist, 

Scheitelpunktbild {V(i,j)}, bei dem alle Randbildpunkte durch einen der folgenden Scheitelpunktarten iden- 
tifiziert sind: isolierter Punkt, Endpunkt, Randpunkt, T-Verbindung, X-Verbindung, Kurvenpunkt oder Ecke, 
Segmentbild {S(i,j)} f bei dem alle aus {V(i.j)} erhaltenen Randbildpunkte vom Scheitelpunkttyp "Rand- 
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punkt" der Mark© des entsprechenden Bildpunktes in (L(i,j)} zugeordnet werden und die anderen Bild- 
punkte null sind,{L(i,j)} 

Ableiten von Bildpunkt n mit identischen Marken aus {S(i,j)}, 

Erzeugen von Liniendatenstrukturen zum Zweck der Darst (lung von Linien, wob i di Datenstrukturen Felder 
zum Darstellen anderor linienspezifischer Daten umiassen wie solche, die von darauffolgenden Verfahrens- 
schritten und Prut ungen erfordert werden, 

Speichern jedes Satzes von Koordinaten von identisch markierten Bildpunkten, die aus {S(i,j)} abgeleitet wor- 
den sind, in einer Liniendatenstruktur; 

f Or jede der Liniendatenstrukturen Schatzen der Lange der entsprechenden Linie, 

- Eliminieren von Liniendatenstrukturen, die alien Linien entsprechen, deren Lange geringer ist als ein vorbe- 
stimmter Wert, aus den Liniendatenstrukturen, 

- Bestimmen des senkrechten Radius und der senkrechten Ausrichtung jeder Linie durch lineare Regression 
und Aufzeichnen dieser Werte in den Liniendatenstrukturen als Attribute der entsprechenden Linie, 
Eliminieren aller Liniendatenstrukturen, f Or die die durch lineare Regression bestimmte Konformitat mit einer 
geraden Linie unterhalb eines gewissen vorbestimmten Wertes liegt, 

lineares Extrapolieren jeder nichteliminierten Linie von jedem Ende aus, bis sie entweder die Bildgrenze Oder 
einen Randpunkt nach Bestimmung durch {S(i,j)}, der nicht zu dieser Linie gehort, schneidet, und Aufzeichnen 
der Positionen beider Schnittpunkte als Attribute der entsprechenden Linie in den Liniendatenstrukturen. 

Verfahren nach Anspruch 3, bei dem die Identifizierung des Scheitelpunkts vom Typ ECKE folgendes umfaBt: 

Krummungsbestimmung von verbundenen Ketten aus Randbildpunkten mit Ausnahme von Bildpunkten, bei 
denen es sich um Verbindungen Oder Endpunkte handelt, 
ortliche nichtmaximale Unterdruckung der Krummungswerte, 

Schwellwertbildung der verbleibenden, von null abweichenden Krummungswerte, wobei lediglich solche Punk- 
te als Ecken angenommen werden, deren Werte Ober einem gewissen Schwellwert liegen. 

Verfahren nach Anspruch 3, bei dem die Identifizierung des Scheitelpunkts vom Typ KURVE folgendes umfaBt: 

Krummungsbestimmung von verbundenen Ketten aus Randbildpunkten mit Ausnahme von Bildpunkten, bei 
denen es sich um Verbindungen Oder Endpunkte handelt, 

Schwellwertbildung der verbleibenden KrOmmungswerte, wobei lediglich solche Punkte als Kurven angenom- 
men werden, deren Werte Ober einem gewissen Schwellwert liegen. 

Verfahren nach Anspruch 3, bei dem die gewichtete lineare Regression folgendes umfaBt: 

Berechnung von gewichteten positionsmaBigen Momenten, bis zur zweiten Ordnung, von Randbildpunkten 
in jeder Liniendatenstruktur, wobei die Gewichtung jedes Bildpunktbeitrages bei (i,j) eine Funktion der Gradi- 
entengroBe G(i,j) ist, 

- Verwendung der gewichteten Momente zur Bestimmung der Parameter der geraden Linie mit der besten An- 
passung zur Minimierung der gewichteten Summe quadrierter Fehler (von Bildpunktpositonen zu der para- 
metrisierten Linie). 

Verfahren nach Anspruch 2, bei dem es sich bei den Grundelementen der Zwischenstufe um Liniengruppen han- 
delt, die durch einen linear angehauften BOndelungsprozeB erhalten worden sind. 

Verfahren nach Anspruch 7, bei dem in dem linear angehauften BOndelungsprozeB zwei beliebige Linien Oder 
zwei Liniengruppen Oder eine Linie und eine Liniengruppe in eine einzige Liniengruppe zusammengef Ohrt werden, 
wenn: 

sie sich in ihrer Ausrichtung um weniger als einen gewissen vorbestimmten Wert unterscheiden, 
der Abstand bei ihrer nachsten Annaherung unter einem gewissen vorbestimmten Wert liegt, 
sie uberlappende Extrapolierungen aufweisen, 

das Ergebnis der F-Statistik auf der Grundlage der Varianzanalyse derart ist, daB die zwei in Betracht gezo- 
genen Bildpunktverteilungen, die am wahrscheinlichsten sind, zu einer einzigen linearen Bildpunktverteilung 
gehoren, die das zusammengef Oh rte Bundel ausmacht. 
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9. V rfahren nach Anspruch 7, bei dem in Frage kommende Hypothesen fur die Signal-Schatten-Grenz aus jeder 
moglichen Liniengrupp , jedem moglich n Paar von Liniengruppen, jedem moglichen Triple von Liniengruppen 
und jedem m~glich n Quartett von Liniengruppen gebildet werden, und zwar auBerdem inschlieBlich d r Hypo- 
these, daB k ine Signal-Schatten-Gr nze voriiegt. 

s 

10. Verfahren nach Anspruch 7, bei dem die VerifikationsprOfungen auf hoher Stufe auf die Annehmbarkeit einer hy- 
pothetischen Signal-Schatten-Grenze al!e beliebigen togischen Kombinationen von Prut ungen auf folgende Punk- 
te umfassen: 

10 - GrenzschlieBung, wobei bei der betrachteten Hypothese, um diese PrOfung zu bestehen, die extrapolierten 

Liniengruppen jeweils an einer anderen Liniengruppe in der gleichen Hypothese Oder am Rand des Bildes 
enden mussen oder die Liniengruppe bzw. den Rand schneiden mussen, 

Grenzform, wobei bei der betrachteten Hypothese. um diese PrOfung zu bestehen, die Grenzform darauf 
eingeschrankt ist, rechtwinklig oder beinahe rechtwinklig zu sein, oder im Fall, daB in der Hypothese weniger 
is als vier Seiten vorliegen, die sichtbaren Seiten eingeschrankt sind, um mit einer teilweise sichtbaren recht- 

winkligen oder beinahe rechtwinkligen Form Obereinzustimmen, 

Signal-Schatten-Kontrast, wobei bei der betrachteten Hypothese, um diese Prufung zu bestehen, der Kontrast 
eingeschrankt ist, so daB er im groBen und ganzen mit der Tatsache ubereinstimmt, daB in dem Bereich von 
{X(i,j)}, der dem Schattenbereich entspricht, ein Grad an Belichtung voriiegt, der niedriger ist als der Belich- 
20 tungsgrad in dem Bereich von {X(i,j)}, der dem Signalbereich entspricht, 

- Grenzort, wobei bei der betrachteten Hypothese, um diese PrOfung zu bestehen, die Mitte des Bildes entweder 
im Signalbereich enthalten sein muB oder nahe an einer Signal-Schatten-Grenze liegen muB, 
Grenzausrichtung, wobei bei der betrachteten Hypothese, um diese Prufung zu bestehen, jede sichtbare Seite 
des Rechtecks parallel zu einer Seite des Bildes verlaufen muB, 

25 - Grenzsymmetrie, wobei bei der betrachteten Hypothese, um diese Prufung zu bestehen, zwei oder vier sicht- 

bare Seiten vorliegen mussen, die ausgerichtet sind (wieoben) undvom Rand des Bildes ungefahraquidistant 
sind, 

- Grenzvollendung, wobei bei der betrachteten Hypothese, um diese PrOfung zu bestehen, der Anteil der Grenz- 
umfangslange, fur den Befund von den Grundelementen niedriger Stufe voriiegt, uber einem gewissen 

30 Schwellwert liegen muB. 

11. Verfahren nach Anspruch 2, wobei jeder Hypothese Kosten zugewiesen werden, wobei diese Kosten am niedrig- 
sten sind, wenn es sich bei der Hypothese um ein Rechteck mit vier sichtbaren Seiten handelt, und wobei die 
Kosten mit der Abnahme der Anzahl der sichtbaren Seiten eines Rechtecks in einer Hypothese ansteigen und die 

35 Kosten am hochsten sind, wenn die Hypothese keine sichtbaren Seiten aufweist. 

12. Verfahren nach Anspruch 11 r zum weiteren Verfeinern der Kosten jeder Hypothese, wobei die Kosten gemaB 
demjenigen Anteil der Signal-Schatten-Grenze, der nicht durch Befund von den Grundelementen niedriger Stufe 
gestutzt wird, ansteigen. 

40 

13. Verfahren nach einem der vorhergehenden AnsprOche, wobei die Bildverarbeitung, die an der digitalen Signaldar- 
stellung des Rontgenbildes durchgefOh rt wird, auf Bildpunktwerte des Bildes innerhalb der Signal-Schatten-Grenze 
begrenzt ist. 

45 14. Bei einem Verfahren zum Lesen eines in einem photostimulierbaren Leuchtstoffschirm gespeicherten Strahlungs- 
bildes durch Abtasten des Schirms mit stimulierender Strahlung, Erfassen von bei der Stimulierung abgehendem 
Licht und Umwandeln des erfaBten Lichts in eine digitale Signaldarstellung, ein Verfahren zum Bestimmen einer 
Signal-Schatten-Grenze gemaB einem der vorhergehenden AnsprOche. 

50 

Revendications 

1. Procede de determination de I'emplacement de la frontiere signal/ombre dans une image radiograph ique repre- 
sentee dans une representation de signal numerique, comprenant les etapes consistant a: 

55 

i) Extraire des primitiv s de bas niveau de I'imag radiograph ique {X(ij)}, 

ii) Former un nombre reduit de primitives de niveau intermediaire a partir des primitives de bas niveau, 

iii) Elaborer des hypotheses comme solutions possibles relatives a I'emplacement de la frontiere signal-ombre, 
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k partir do combinaisons de primitives de niveau intermediate, etape durant laquelle chaque combinaison st 
soumise a des tests de niveau intermediate, k (a suite desqu Is des hypotheses parttelles ou completes sont 
rejete sou accepts s, 

iv) Ftealiser des tests de verification de haut niveau sur chaque hypothese, a ia suite desquels des hypotheses 
sont rejet6es ou acceptees k un certain coGt, ledit cout d6finissant un classement de ['acceptability d'une 
hypothese, 

v) S6lectionner I'hypothese de moindre coOt. 

Proc6d6 selon la revendicatbn 1, dans lequel lesdites primitives de bas niveau sont des lignes. 

Procede selon la revendicatbn 2, dans lequel lesdites tignes sont extraites par realisation des Stapes suivantes 
sur \a representation de signal num6rique: 

sous-6chantiltonnage, 

- realisation d'une detection de bords sur ladite image sous-echanti!lonn6e pour produire les images intrinsfc- 
ques suivantes: 

- carre de I'amplitude de gradient local {G(i,j)J, 

- orientation locale de bord (0(i,j)}, 

image de bord {E(i,j)J, dans laquelle tous les pixels dont le carre de I'amplitude de gradient est inferieur 
k une valeur seuii pr6d6terminee sont fixes k zero, et dans laquelle tous les pixels de bord qui n'atfectent 
pas la connectivite locale d'autres pixels de bord ont ete supprimes par une operation d'allegement, 
image de designation {L(i, j)}, dans laquelle tous les pixels de bord mutueltement connects de I'image de 
bord aliegee partagent une designation commune, quoique autrement unique, 

image d'orientatbn locale de bord, obtenue par rdgressbn lineaire de pixels de bord connectes dans un 
voisinage local, 

image de sommets {V(i,j)}, dans laquelle tous les pixels de bord sont identifies comme etant Tun des types 
de sommet suivants: point isoie, terminaison, point de bord, jonction en T, jonction en X, point de courbe, 
ou coin, 

image de segment {S(i,j)}, dans laquelle on attribue k tous les pixels de bord du type de sommet 'point 
de bord', tels qu'obtenus a partir de {V(i,j)}, la designation du pixet correspondant dans {L(i,j)J, et les autres 
pixels ont pour valeur zero, 

extraction de {S(i, j)} des pixels ayant des designations identiques, 

creation de structures de donn6es de lignes dans le but de repr6senter des lignes, lesdites structures de 
donn6es comprenant des champs pour representor d'autres donnees specifiques k des lignes, necessaires 
pour des etapes et des tests ult6rieurs du procede, 

memorisation de chaque ensemble de coordonnees de pixels d6signes de maniere identique, extraits de S{ 
(i,j)}, dans une structure de donnees de lignes, 

estimation, pour chacune des structures de donnees de lignes, de la longueur de la ligne correspondant e, 
elimination desdites structures de donn6es de lignes les structures de donnees de lignes correspondant a 
toutes les lignes dont la longueur est interieure k une certaine valeur predetermin6e, 

- determination du rayon normal et de i'orientation de chaque ligne par regression lineaire et enregistrement 
de ces vaieurs dans lesdites structures de donnees de lignes en tant qu'attributs de la ligne correspondante, 
elimination de toutes les structures de donnees de lignes pour lesquelles la conform ite par rapport a une 
drorte, telle que d6termin6e par regression lineaire, est interieure a une certaine valeur predetermine, 
extrapolation lineaire de chaque ligne non 6!imin6e depuis chaque extremite, jusqu'a ce qu'elle coupe soit la 
frontiers de I'image, soit un point de bord tel que determine par {S(i,j)} ( qui n'appartient pas k cette ligne, et 
enregistrement, dans les structures de donn6es de lignes, des positions des deux points d'intersection en tant 
qu'attributs de la ligne correspondante. 

Procede selon la revendicatbn 3, dans lequel ladite identification du type de sommet COIN comprend 

la determination de la courbure de chaines connectees de pixels de bord, k I'exclusion des pixels qui sont 

des jonctions ou des pixels terminaux, 

la suppresion locale non maximale des vaieurs de courbure, 

la fixation d'un seuil pour les vaieurs de courbure non nulles restantes, en acceptant en tant que coins uni- 
quement les points dont la valeur est sup6rieure k un certain seuil. 
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5. Proc6d6 seton la revendication 3, dans lequel tadite identification du typ de sommet COURBE comprend 

la determination d la courbure de chatnes connectees de pixels de bord, a I'exclusion d s pix Is qui sont 
des jonctions ou des pixels terminaux, 
s la fixation d'un seuil pour les valeurs de courbure, en acceptant en tant que courbes uniquement les points 

dont la valeur est superieure a un certain seuil. 

6. Proc6d6 selon la revendication 3, dans lequel ia regression Iin6aire pond6r6e comprend: 

10 - |e calcul de moments de position pond6r6s, jusqu'au deuxieme ordre, de pixels de bord dans chaque structure 

de donn6es de ligne, dans lequel la ponderation de chaque contribution de pixel en (i.j) est une fonction de 
Pamplitude de gradient G(i,j), 

I'utilisation des moments pond6r6s pour determiner les parametres de la droite de meilleur ajustement de 
maniere a minimiser la somme pond6ree du carre des erreurs (des positions des pixels par rapport a la ligne 
is parametree). 

7. Proc6d6 selon la revendication 2, dans lequel lesdites primitives de niveau intermediate sont des groupes de 
lignes obtenus par une methode de mise en faisceaux Y\n6aire agglom6rante. 

so 8. Procede selon la revendication 7, dans lequel, dans la methode de mise en faisceaux lineaire agglomerante, deux 
lignes ou deux groupes de lignes, ou une ligne et un groupe de lignes quelconques sont fusionnes en un groupe 
de lignes unique si: 

leur difference d'orientation est inferieure a une certaine valeur predetermines, 
leur distance d'approche minimale est interieure a une certaine valeur pred6termin6e, 
ils ont des extrapolations en chevauchement 

le resultat de la statistique F base sur I'anafyse de la variance est tel que les deux distributions de pixels 
consid6r6es appartiennent tres probablement a une distribution de pixels lineaire unique qui constitue le fais- 
ceau fusionne. 

Proc6de selon la revendicatbn 7, ou des hypotheses candidates pour la frontiere signal/ombre sont 6labor6es a 
partir de chaque groupe de lignes, paire de groupes de lignes, triplet de groupes de lignes, et quadruplet de 
groupes de lignes possibles, y compris, en outre, i'hypothese qu'il n"y a aucune frontiere signal/ombre. 

35 10. Proc6de selon la revendication 7, dans lequel lesdits tests de verification de haut niveau d'acceptabilite d'une 
frontiere signal/ombre hypothetique, comprennent toute combinaison logique des tests suivants de: 

- f ermeture de la frontiere; pour qu'un tel test sort passe avec succes, pour I'hypothdse considerde, fes groupes 
de lignes extrapoies doivent chacun aboutir sur, ou couper, un autre groupe de lignes dans la m§me hy pothese, 

40 ou le bord de I'image, 

forme de la frontiere; pour qu'un tel test soit pass6 avec succ&s, pour I'hypothese consid6r6e, on contraint ia 
forme a etre rectangulaire, ou pratiquement rectangulaire ou, dans le cas ou il y a moins de quatre c6t6s dans 
I'hypothese, on contraint les cot6s visibles a §tre compatibles avec une forme rectangulaire ou pratiquement 
rectangulaire partiellement visible, 

45 - contraste signal/ombre; pour qu'un tel test soit passe avec succes, dans Thypothese considered on contraint 

le contraste k etre largement compatible avec le fait qu'il y ait un niveau d'exposition dans la region de {X(i, 
j)} correspondant a la region d'ombre, inferieur au niveau d'exposition dans la region de {X(i,j)} correspondant 
a la region de signal, 

emplacement de la frontiere; pour qu'un tel test soit pass6 avec succ6s, pour I'hypothese consid6r6e, le centre 
so de I'image doit etre soit contenu k I'int6rieur de la region de signal, soit §tre proche d'une frontiere signal/ombre, 

alignement de la frontiere; pour qu'un tel test soit passe avec succes, pour I'hypothese consider, chaque 
cote visible du rectangle doit etre paralieie a un cote de Pimage, 

symetrie de la frontiere; pour qu'un tel test soit passe avec succes, pour I'hypothese consider6e, il doit y avoir 
deux ou quatre c6\6s visibles, qui soient alignes (comme ci-dessus), et approximativement equidistants du 
55 bord de I'image, 

achevement de la frontiere; pour qu'un tel test soit passe avec succes, pour I'hypothese consider6e, la fraction 
de la longueur du perimetre de la frontiere, pour laquelle il existe des preuves provenant des primitives de 
bas niveau, doit etre superieure a un certain seuil. 
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11. Precede s Ion la revendication 2, dans lequ I un coOt est attribue a chaque hypothes , lequel coOt est te moins 
6\ev6 si I'hypothese est celle d'un rectangle ayant quatre cdtes visibles, et lequel coOt augment a mesure que le 
nombr de c6tes visibl s d'un rectangle dans un hypothese diminue, 1 le coOt est le plus eleve si I'hypothese 
n'a aucun cdt6 visibl . 

12. Precede selon la revendication 11 , en vue d'affiner le coOt de chaque hypothese, dans lequel le cout est augments 
seloh la traction de la frontiere signal/ombre qui n'est pas appuyee par des preuves provenant des primitives de 
bas niveau. 

13. Precede selon Tune quelconque des revendtcations precedentes, dans lequel te traitement de ('image realise sur 
ladite representation de signal numenque de ladite image radiograph ique est limits aux valeurs de pixels de ladite 
image a I'interieur de la frontiere signal/ombre. 

14. Dans un precede de lecture d'une image de rayonnement memorised dans un ecran luminescent photostimulable 
par balayage dudit ecran avec un rayonnement de stimulation, detection de la lumiere 6mise suite a la stimulation, 
et conversion de la lumiere d6tectee en une representation d'image numerique, precede de determination d'une 
frontiere signal/ombre selon I'une quelconque des revendications prec6dentes. 
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FIG. 9 
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